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SUMMARY/OVERVIEW:  Distribution  Unlimited 

Studies  were  conducted  in  several  relevant  areas,  including  (1)  validation  of  the  chemistry  and 
transport  models  against  the  extinction  of  ultra-lean  premixed  hydrogen-air  mixmres,  (2)  a 
comprehensive  theoretical  analysis  of  the  reaction  kinetics  of  carbon  monoxide  and  the  hydroxyl 
radical,  (3)  a  theoretical  kinetic  study  of  the  decomposition  of  ethylene  oxide;  (4)  a  gas-kinetic 
analysis  for  the  transport  properties  of  long  chain  molecules  in  dilute  gases,  (5)  quantum-chemistry, 
master  equation  modeling  of  the  unimolecular  decomposition  of  ortho-benzyne,  (6)  extension  of  the 
previously  developed  hydrogen/ carbon  model  to  combustion  pressures  as  high  as  600  atm,  (7)  an 
updated  kinetic  mechanism  of  small-hydrocarbon  fuel  combustion  for  use  as  a  kinetic  foundation  of 
higher  hydrocarbon  combustion,  and  (8)  a  methodology  for  kinetic  uncertainty  propagation.  These 
projects  represent  the  two  key  ingredients  for  meeting  the  overall  project  objectives:  (a)  an  accurate 
physico-chemical  property  database  for  combustion  kinetics,  and  (b)  a  unified  and  optimized  kinetic 
model  for  liquid  aliphatic  and  aromatic  fuel  combustion  with  quantifiable  uncertainties. 

TECHNICAL  DISCUSSION: 

1.  Validation  of  Chemistry  and  Transport  Models  for  H^/CO  Combustion 

In  this  collaborative  project  (with  Professor  Fokion  Egolfopoulos  of  USC),  we  validated  a  Hj- 
CO  kinetic  model  and  its  associated,  revised  transport  database  [1]  against  the  extinction  of 
premixed  ultralean  H2/air  flames  [2].  Theoretical  attention  has  been  placed  on  the  relative  sensitivity 
of  predicted  flame  extinction  with  respect  to  reaction  kinetics  and  species  diffusion  coefficients. 

Under  the  conditions  considered  in  this  smdy,  the  sensitivity  of  diffusion  coefficients  on  the 
extinction  response  was  found  to  be  significant  and  for  certain  species  this  sensitivity  is  greater  than 
the  kinetics  rate  constants.  This  study  demonstrates  the  importance  of  an  accurate  transport  library 
in  the  development  of  predictive  reaction  kinetic  models. 

This  work  is  reported  in  Ref  2. 
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2.  Theoretical  Analysis  of  the  CO  +  OH  — »  CO^  +  H  Reaction 

In  a  previous  study  [1],  we  demonstrated  that  the  currendy  available  combustion  data  of  and 
CO  combustion  cannot  be  reconciled  unless  a  bi-exponential  rate  expression  is  used  for  the  reaction 
CO  +  OH  — >  CO2  +  H.  While  experimental  evidence  supports  this  bi-exponential  temperature 
dependence,  only  a  limit  number  of  kinetic  models  employ  such  a  rate  expression.  Theoretically  the 
characteristic  dependence  of  the  rate  constant  on  temperature  as  well  as  the  pressure  dependence 
remains  to  be  carefully  examined.  In  the  present  study,  the  rate  coefficients  of  CO  +  OH  -> 
products,  including  the  chemically  activated  addition/dissociation  channel  to  produce  COj  +  H  and 
the  addition/stabilization  reaction,  are  analyzed  with  RRKM/master  equation  analyses.  The  master 
equation  of  collision  energy  transfer  was  solved  with  a  newly  implemented  Monte  Carlo  code.  The 
RRKM  analyses  are  based  on  the  CCSD(T)/cc-pvTZ  potential  energy  surface  of  Yu  et  al.  (Chem. 
Phys.  Lett.  2001,  349,  547-554).  It  is  shown  that  over  the  temperamre  range  of  80  to  2500  K  and 
pressures  from  1  Torr  to  800  Bar  the  experimental  data  can  be  satisfactorily  reproduced  by  lowering 
the  CCSD(T)/cc-pvTZ  energy  barrier  for  the  CO2  +  H  exit  channel  by  1  kcal/mol  and  more 
importantly,  by  considering  an  equilibrium  factor  in  the  thermal  rate  constant  formulation.  This 
factor  accounts  for  the  populations  of  rovibrationally  excited  trans-  and  ar-HOCO,  which  are 
allowed  to  dissociate  only  through  specific  paths  that  are  open  to  them.  Figure  1  shows  the  rate 
constant  of  CO  +  OH  — »  products  for  selected  temperatures  over  a  wide  range  of  temperature.  It 
is  seen  that  the  current  RRKM/master  equation  model  predicts  very  well  the  pressure  dependent 
rate  constant  over  a  wide  range  of  temperatures. 

By  modeling  the  isothermal  but  pressure-dependent  rate  data  of  FuUe  et  al.  (J  Chem  Phys  1996, 
105,  983-1000)  over  the  temperature  range  from  98  to  819  K,  we  obtained  an  value  equal  to 

150  cm"’  for  M  =  He.  Figure  2  shows  the  comparison  of  these  pressure  dependent  data  with  the 
predicted  rate  constants.  The  values  for  M  =  Nj,  Ar,  CF4  and  SF^  were  also  obtained  by 

fitting  the  OH  and  OD  data  at  298  K.  Based  on  the  current  theoretical  analyses,  we  recommended 
that  the  following  low-pressure  limit  rate  expression  be  used  for  CO  +  OH  — >  CO2  +  H, 

0  ^cm^molecule~'s~' j  =  1.17xlO~^^T^'^^^  exp^-^^  +  9.56xl0’*^^T'*°  ^^‘*  exp  — —  . 

The  above  expression  is  valid  for  temperature  from  120-2500  K.  In  addition,  for  P  (Bar)  < 
9xl0''^^T^'^  exp(520/T) ,  the  rate  constant  is  practically  at  the  low  pressure  limit.  FaU-off 
parameterization  is  also  proposed  over  the  entire  range  of  pressure  for  the  rate  coefficients  of  CO  + 
OH  ->  CO2  +  H  and  CO  +  OH  —>  HOCO  over  the  temperature  range  from  120  to  2500  K.. 

For  800  K  <  T  <  2500  K  the  theoretical  rate  expression  given  above  agrees  remarkably  well  with 
the  same  rate  obtained  by  optimizing  our  recent  H2/ CO  model  [1]  against  the  available  combustion 
data  (<5%  difference).  In  addition,  the  current  theoretical  analysis  shows  that  for  combustion 
applications  the  pressure  dependence  of  the  rate  constant  for  CO  +  OH  —>  CO2  +  H  is  negligible 
since  for  T  >  1200  K  the  low-pressure  expression  given  above  is  valid  for  pressure  <  200  Bar. 

This  work  is  reported  in  Ref  3. 
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Figure  1.  Calculated  M  =  He)  and  obsen'cd  (other  symbols)  pressure  dependence  of 
the  rate  constant  of  CO  +  OH  — ♦  produce  at  (a)  190  K  (V  ;  M  =  He,  T=:210K;^:M  = 
He,  r=  220  K;  +  :  M  =  At,  T=  216  K;  ♦:  M  =  He),  (b)  250  K  (V ;  M  =  He,  interpolated 
to  T=  250  K;  A:  M  =  He,  T  =  240  KiX:  M  =  At,  T=  251  K;^:  M  =  He),  (c)  298  K  for 
M  =  He,  (d)  400  K  for  M  =  He  (P:  T=  396  K;  V :  interpolated  to  T=  400  K;  A :  T=  373 
K;  ■:  interpolated  to  T  =  400  K),  (b)  500  K  (A :  M  =  He,  T=  495  K;  D:  M  =  He,  T=  523 
PC;  V ;  M  =  He,  extrapolated  to  T  =  500K;  ■ :  M  —  He,  extrapolated  to  T  =  500K;  ^ :  M  = 
Ar.  T  =  511  K;  •:  M  =  Ar,  T  =  511  K;  ♦ :  M  =  He),  and  (c)  819  K  P:  M  =  He, 
interpolated  to  T  =  819  K;  :  M  =  Ar,  T  =  800  K;  ® :  M  =  Ar,  interpolated  to  T  =  819  K ; 
♦:  M  =  He1. 
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Figure  2.  Calculated  (^ ,  M  =  He)  and  observed  (other  symbols) 
temperature  and  pressure  dependence  of  CO  +  OH  — »  products.  D;  M 
=  Ar  at  ~2  Bar  +;  M  =  Ar  at  0.19  -  0.83  Bar;  =  Ar  at  50  Torr; 

M  =  Ar  at  1  Bar;  A;  M  =  He  at  1±0.2  Bar;^:  M  =  He  at  10±1 
Bar;®:  M  =  He  at  80±15  Bar;^  :  M  =  He  at  650  Bar;  vibration 
relaxation  rate  constant  of  OH(l>  ~1)  +  CO  •->  OH(l>=2)  +  CO 
(approximately  equal  to  >fc«). 


3.  Reaction  Kinetics  of  Ethylene  Oxide  Decomposition 

The  unimolecular  decomposition  of  ethylene  oxide  and  the  oxiranyl  radial  is  examined  by 
molecular  orbital  calculations,  Rice-Ramsperger-Kassel-Marcus  (RRKM) /Master  Equation  analysis, 
and  detailed  kinetic  modeling  of  ethylene  oxide  pyrolysis  in  a  single-pulse  shock  tube.  The  study 
was  motivated  by  the  need  to  better  understand  the  kinetics  of  reactions  involving  spin  state 
crossing.  The  development  of  this  model  was  also  made  in  response  to  the  need  of  the  Wright- 
Patterson  research  lab  for  a  reliable  reaction  model  for  ethylene  oxide  combustion. 

It  was  found  that  the  largest  energy  barrier  to  the  decomposition  of  ethylene  oxide  lies  in  its 
initial  isomerization  to  form  acetaldehyde,  and  in  agreement  with  previous  studies,  the  isomerization 
was  found  to  proceed  through  the  •CH2CH2O  •  biradical.  Because  of  the  biradical  nature  of  the 
transition  states  and  intermediate,  the  energy  barriers  for  the  initial  C-O  rupture  in  ethylene  oxide 
and  the  subsequent  1,2-H  shift  remain  highly  uncertainty  (see.  Figure  3).  On  the  other  hand,  an 
overall  isomerization  energy  barrier  of  59  ±  2  kcal/mol  was  found  to  satisfactorily  explain  the 
available  single  pulse  shock  tube  data  (Figure  4).  This  barrier  height  appears  to  be  in  line  with  the 
estimates  made  from  an  approximate  spin-corrected  procedure  at  the  MP4/6-31-iG(d)  and 
QCISD(T)/6-31G(d)  levels  of  theory. 

The  high  pressure  limit  rate  constant  for  the  unimolecular  decomposition  of  ethylene  oxide  was 
calculated  as 

(l/s)  =  3.74  X 10’"  for  600  <  T  <  2000  K. 


3 


The  dominant  channel  for  the  unimolecular  decomposition  of  ethylene  oxide  was  found  to  form 
CHj  +  HCO,  and  at  the  ambient  pressure  and  for  T  >  800  K  it  accounts  for  >90%  of  the  total  rate 
constant. 

The  present  work  constitutes  a  step  forward  in  our  understanding  of  the  singlet  manifold  of  the 
C2H4O  potential  energy  surface.  A  combination  of  this  singlet  PES  with  that  of  the  triplet  states  is 
required  to  address  the  chemistry  of  the  O  +  C2H4  reaction.  This  combination  of  the  surfaces 
requires  calculations  of  the  spin-orbit  coupling  between  the  singlet  and  triplet  ring-opened 
•CH2CH2O  •  biradicals.  We  are  hopeful  that,  with  the  singlet  surface  now  well  understood,  we  will 


be  able  to  accurately  predict  the  branching  rations  and  the  respective  rate  constants  for  the  reaction 
of  O  +  C2H4. 

This  work  is  reported  in  Ref  4. 

4.  Transport  Theory  of  Long  Chain  Molecules  in  Dilute  Gases 

One  of  the  major  theoretical  hurdles  to  obtain  a  quantitative  understanding  of  the  combustion 
chemistry  of  jet  fuels  and  their  surrogates  is  the  lack  of  a  suitable  theory  to  describe  the  diffusion 
coefficient  of  long  chain  aliphatic  molecules.  Although  the  Chapman-Enskog  theory  is  traditionally 
employed  to  calculate  the  binary  diffusion  coefficient,  its  underlying  assumption  of  spherical 
potentials  for  the  collision  pair  renders  the  applicability  of  the  theory  questionable.  In  this  work,  a 
rigorous,  gas-kinetic  theory  analysis  was  carried  out,  and  analytical  equations  are  developed  for  the 
diffusion  coefficient  of  long-chain  molecules  in  the  limits  of  specular  and  diffuse  scattering. 

A  long-chain  molecule  may  be  approximated  as  a  straight  and  rigid  cylinder.  Consider  a  cylinder 
with  length  equal  to  E  and  diameter  D.  Gas-kinetic  theory  analysis  shows  that  the  instantaneous 
drag  force  due  to  relative  motion  of  a  cylinder  in  a  dilute  gas  may  be  given  by 

F,„  =  1^(1  +  5  sin^  a/4)  2^«!,kTNLDVQ!-y]-  (1) 
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Figure  4.  Compari.con  of  predicted  (lines)  and  experimental  (symbols,  [8]) 
concentration  profile  of  ethylene  during  shock  tube  pyrolysis  of  ethylene  oxide 
(dwell  time  =  2  ms).  Computations  were  carried  out  for  Mixture  A  (solid  lines) 
and  Mixture  C  (dashed  line.s). 


=  61  kcal/mol 


A:  1%oCjH,0-99%Ar,  C5  =  2.3><10'  mol/cirf 
B:  0.25%c-CjH,O-99.75»/nAr,  C5  =  5.2x10“  mol/cm“ 
C:  0.25%c.C2H^O-99.75%Ar,  C5  =  2.2x10“  mol/cirf 
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where  the  subscript  denotes  specular  or  diffuse  scattering,  a  is  the  angle  of  the  cylinder 

relative  to  the  flow,  is  the  reduced  mass,  k  is  the  Boltzmann  constant,  T  is  the  temperature,  N  is 
the  gas  number  density,  V  is  the  relative  velocity.  In  eq.  (1)  is  the  reduced  collision  integral 

given  by 


no 

LD 


(2) 


where  is  the  collision  cross  section,  defined  in  the  two  limiting  scattering  scenario  as 

=  ^  r  0-cos;ir)^^, 

•1-00 

and 


(3a) 


db 


1  +  i  ^3in^ 

+  ^^{\-coiX)dh+  £“(1  -  cos  x) 


(3b) 


respectively.  In  eqs  (3),  the  angle  X  is  given  by 

X{i 


,b)  =  n-2b  P - . 

J.  ,3 


(4) 


where  ris  the  center-to-axial  distance  between  the  gas  molecule  and  cylinder,  is  the  distance  of  the 
closest  encounter,  and  0(r,L)  is  the  potential  function  between  the  gas  molecules  and  cylinder.  In 
the  limit  of  rigid-body  coUision,  it  may  be  shown  that 

^  4/3 

and 

Qti’')*  =l-t-3;rV64.  (5b) 

Realistically  molecular  scattering  always  lies  between  the  two  limiting  scattering  models.  Similar  to 
the  treatment  of  nanoparticle  transport,  a  parameter  may  be  introduced  here  to  give  a  mixed 
scattering  model,  such  that  an  average,  reduced  collision  integral  may  be  defined  as 


where  ^  is  termed  as  the  momentum  accommodation  function,  which  ranges  from  0  to  1.  From 
the  knowledge  of  gas-nanoparticle  scattering,  it  may  be  speculated  that  for  long-chain  aliphatic 
molecules  (p  is  close  to  unity  [5]. 

The  instantaneous  binary  diffusion  coefficient  of  a  cylinder  in  a  dilute  gas  may  be  obtained  from 
the  Stokes-Einstein  relation,  through  the  drag  coefficient,  as 
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(7) 


A2=2 


kT 


(l  +  5sin^ 


An  angle-averaged  diffusion  coefficient  is  obtained  by  replacing  a  by  the  mean  angle  OC  of  the 
orientation  of  the  cylinder  with  respect  to  its  diffusion  direction. 

The  remaining  work  will  be  (a)  to  examine  the  potential  function  0(r,L)  and  to  determine  a 
suitable,  parameterized  form  for  this  function,  and  (b)  to  carry  out  suitable  experiments  to  determine 
the  diffusion  coefficients  of  representative,  long-chain  ahphatic  molecules.  Both  components  of  the 
study  are  required  to  validate  the  theory  developed. 


5.  Kinetics  of  unimolecular  decomposition  of  or^'^o-benzyne 

The  kinetics  of  o^Ao-benzyne  is  an  integral  part  of  aromatics  oxidation  at  high  temperatures. 
Previously,  it  was  concluded  that  the  dominant  channel  of  ort/)o-henzyne  decomposition  is  through 
the  concerted  unimolecular  dissociation  leading  to  diacetylene  and  acetylene, 

+  (8) 

Recent  shock-tube  studies  at  DLR  suggest  that  an  additional  reaction  channel  may  become  active  at 
temperatures  >  1600  K.  This  channel  is  characterized  by  the  production  of  the  H  atom, 

o-QH,-*QH,  +  2H  (9) 

The  shock  tube  study  used  1,2-diiodobenzene  is  the  precursor  for  or^/>o-benzyne. 

In  this  study,  we  examined  the  potential  energy  surface  of  orf/)o-henzyne  decomposition  by 
molecular  orbital  calculations.  The  computation  was  carried  out  at  the  CCSD('I)/6- 
31G(d)//B3LYP/6-31G(d)  level  of  theory.  The  minimum  energy  barrier  for  H  production  was 
computed  to  be  106  kcal/mol  from  orf/>o-henzyne.  This  energy  barrier  renders  direct  H  production 
implausible  below  2000  K.  Subsequendy,  the  reaction  rate  coefficient  was  analyzed  by  master 
equation  modeling.  Kinetic  modeling  shows  that  below  2000  K  the  apparent  production  of  the  H 
atom  is  due  to  the  decomposition  of  the  precursor  compound,  rather  than  the  decomposition  of 
orth-benzyne.  The  work  is  reported  in  Ref  6. 

6.  HJ CO  combustion  up  to  600  atm 

This  is  a  collaborative  research  project  with  Professor  Brezinsky  at  UIC.  In  this  study,  the  high  • 
pressure  oxidation  of  dilute  CO  mixtures  doped  with  150-200  ppm  of  H2  has  been  studied  behind 
reflected  shock  waves  in  the  UIC  high  pressure  single  pulse  shock  tube.  The  experiments  were 
performed  over  the  temperature  range  from  1000-1500  K  and  pressures  spanning  21-500  bars  for 
stoichiometric  (0=1)  and  fuel  lean  (0=0.5)  oxidation.  Stable  species  sampled  from  the  shock  tube 
were  analyzed  by  standard  GC,  GC/MS  techniques.  The  experimental  data  obtained  were  simulated 
using  a  previously  published  detailed  model  for  H2/ CO  combustion.  Simulations  showed  that  within 
experimental  error  the  kinetic  model  was  able  to  capture  the  experimental  trends  for  the  lower 
pressure  data  sets  (average  nominal  pressures  of  24  and  43  bars).  However  the  model  under  predicts 
the  CO  and  O2  decay  and  subsequent  CO2  formation  for  the  higher  pressure  data  sets  (average 
nominal  pressures  of  256  and  450  bars).  The  current  elevated  pressure  data  sets  span  a  previously 
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unmapped  regime  and  have  served  to  probe  HO2  radical  reactions  which  appear  to  be  among  the 
most  sensitive  reactions  in  the  model  under  these  conditions.  With  updated  rate  parameters  for 
HO2  +  OH  =02  +  H2O  the  model  is  able  to  reconcile  the  elevated  pressure  data  sets  thereby 
extending  the  capability  of  the  reaction  model  to  an  extreme  range  of  conditions. 

This  work  is  reported  in  Ref.  [7]. 


7.  Updated  Reaction  Model 

The  need  to  re-evaluate  our  previous  CJ-C4  kinetic  model  stems  from  the  fact  that  its  H2/ CO 
submodel  had  to  be  revised  because  of  the  recent  downward  revision  of  the  heat  of  formation  of  the 
OH  radical  and  new  kinetic  data  of  H  +  O2  (+M)  — >  HO2  (+M)  made  available  only  recently  [1]. 
Currently  we  have  completed  a  re-evaluation  of  relevant  kinetic  rate  constants  and  thermochemical 
data.  Emphasis  is  placed  on  quantify  the  uncertainty  factors  for  all  rate  coefficients.  Concurrently, 
we  are  exploring  new,  modularixed  optimization  methods,  which  will  allow  future  revision  of  the 
kinetic  submodels  be  made  and  its  optimization  done  in  an  isolated  fashion.  The  resulting  model 
will  have  the  capability  to  examine  uncertainty  propagation,  a  component  critically  needed  to  extend 
the  kinetic  foundation  model  to  large  hydrocarbons. 


8.  Methodology  of  kinetic  uncertainty  propagation 
We  begin  to  develop  a  mathematical  procedure  to 
propagate  kinetic  uncertainty  into  combusting  flow 
simulations.  Figure  1  shows  the  statistical  distributions 
of  CO2  mass  fraction  of  ethylene  oxidation  in  air  in  a 
perfectly  stirred  reactor  (^  =  1.0,  T  =  1400  K  and  p  = 
1  atm)  as  a  function  of  the  residence  time  X.  These 
distributions  were  obtained  by  Monte  Carlo 
simulations.  Although  they  nicely  predict  the  wide 
distributions  under  the  conditions  where  the  residence 
time  is  limited,  and  the  bifurcation  of  the  distribution 
near  the  turning  points,  Monte  Carlo  simulations  are 
inherently  expansive.  Currently,  we  are  looking  at  a 
suitable  mathematical  procedure  by  which  these 
distribution  functions  can  be  computed  directly.  A 
part  of  this  work  is  to  be  reported  in  Ref  [8]. 

Figure  5.  Statistical  distributions  of  PSR  solutions  for  the  CO2 
mass  fraction  (C2H4-air  mixture,  ^  =  1.0,  T  =  1400  K  and  =  1 
atm). 
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EXTRAMURAL-INTRAMURAL  INTERACTIONS 

Our  AFRL  contact  person  is  Dr.  J.  T  Edwards.  Several  visits  to  the  Wright-Patterson  AFRL  were 
made  during  the  project  period.  During  the  visit,  the  PI  discussed  potential  research  collaborations 
with  various  lab  research  personnel,  and  gave  a  seminar  entitled  “Incipient  Soot  Formation  and 
Growth.”  We  also  worked  closely  with  Drs.  Dean  Eklund  and  Skip  Williams  and  delivered  a 
pseudo-detailed  ethylene  oxidation  mechanism  for  Scamjet  simulations. 
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